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Despite more than two decades of intensive in- 
vestigations, the true nature of high tempera- 
ture (high- Tc) superconductivity observed in the 
cuprates remains elusive to the researchers. In 
particular, in the so-called 'underdoped' region, 
the overall behavior of superconductivity devi- 
ates qualitativel'fi^^^ from the standard theoret- 
ical description pioneered by Bardeen, Cooper 
and Schrieffer (BCS).^"^ Recently, the importance 
of phase fluct uatio n of the superconducting or- 
der parameter ^^E^Hm has gained signiflcant sup- 
port from various experimentsJ^^"^ However, the 
microscopic mechanism responsible for the sur- 
prisingly soft phase remains one of the most im- 
portant unsolved puzzles. Here, opposite to the 
standard BCS starting point, we propose a sim- 
ple, solvable low-energy model in the strong cou- 
pling limit, which maps the superconductivity lit- 
erally into a well-understood physics of super- 
fluid in a special dilute bosonic system of local 
pairs of doped holes. In the prototypical material 
(Lai_(5Sr5)2Cu04, without use of any free param- 
eter, a d-wave superconductivity is obtained for 
doping above ~ 5.2%, below which unexpected 
incoherent p-wave pairs dominate. Throughout 
the whole underdoped region, very soft phases 
are found to originate from enormous mass en- 
hancement of the pairs. Furthermore, a striking 
mass divergence is predicted that dictates the oc- 
currence of the observed quantum critical point. 
Our model produces properties of the superfluid 
in good agreement with the experiments, and 
provides new insights into several current puz- 
zles. Owing to its simplicity, this model ofl*ers a 
paradigm of great value in answering the long- 
standing challenges in underdoped cuprates. 

Introduction 

Considering the enormous amount of research activi- 
ties devoted to the problem of high- Tc superconductivity, 
it is hardly an exaggeration to regard it as one of to- 
day's most important unsolved problems in physics. In 
particular, in the underdoped region where only small 
number of 'hole' carriers exist in the system, both the 
low-temperature superconducting phase and the 'nor- 
mal state' above the transition temperature Tc demon- 
strate exotic properties^^^ inconceivable from the stan- 
dard BCS theory. Even if one focuses only on the prop- 
erties of the superfluid, as in this study, still lacking is a 
comprehensive and coherent description capable to simul- 
taneously account for, for example, its d-wave symmetry, 
its rapid thermal depletion, its very soft phase, and its 



dramatic destruction near the 5.2% doping. 

Perhaps the most exotic property of underdoped 
cuprates is the observation of strong diamagnetism^H^, 
one of the signatures of superconductivity, at tempera- 
ture much higher than Tc. This important observation 
suggests the existence of "local" super-current capable of 
screening the magnetic field in the non-superconducting 
state, and th us has been taken as direct evidence of 'pre- 
formed pairs *^^^^^. Following the pioneering observa- 
tion of proportionality between Tc and zero-temperature 
superfiuid density^, ns{T = 0), Emery and KivelsorP 
made a general macroscopic argument: At low carrier 
density the phase of the complex superconducting order 
parameter is soft enough such that its fiuctuation can 
destroy superconductivity before its magnitude vanishes. 
The observed Tc is thus no longer controlled by the pair- 
ing strength U as in any BCS-like mean-field theories, 
but instead by the phase 'stiffness'. 

This important realization, however, left behind the 
most essential microscopic question: "Why is the phase 
in cuprates so soft?" As demonstrated in Fig. la for the 
prototypical high- Tc cuprates (Lai_5Sr5)2Cu04 (LSCO), 
the phase stiffness Vq = h^ns{^)a/Am* ^ evaluated using 
realistic c-axis lattice constant a = 6.6A, effective mass 
771* from twice the experimental mass of the holes, m^ ~ 
4.8me ,^^ ^^ and ns(0) = 5/0.^^ doping level 5 per unit 
cell volume f^o, leads to a T6i ^ 2.2Vb^ much higher than 
the experimental observation. Obviously, a microscopic 
explanation to the very soft phase in the underdoped 
region is one essential piece of the puzzle. 

A related outstanding issue is the nature of the quan- 
tum critical point (QCP) dX 5 = Sqcp ~ 5.2%, the point 
where superconducting phase vanishes at zero tempera- 
ture. As opposed to the above estimation of Tg oc 6, the 
actual Tc diminishes rapidly to zero at Sqcp without an 
immediate appearance of a competing phase. A simple 
picture is still unavailable on exactly what happens to 
the superfiuid near the QCP that suppresses Tc so dra- 
matically, and why the superconductivity ceases to exist 
below this point even at zero temperature. 

Here, we propose a simple low-energy model to ad- 
dress these crucial questions in the underdoped region, 
by starting from the strong coupling limit opposite to the 
BCS approach. The superconductivity is then mapped 
literally into a well-understood problem of superfiuid of 
a special dilute bosonic gas of local pairs of doped holes. 
The model is found to lead to the above puzzling proper- 
ties of the prototypical material (Lai_(5Sr5)2Cu04 with- 
out use of any free parameter. Specifically, a strong com- 
petition with incoherent p-wave symmetry is found to 
enormously enhance the mass of the d-wave pairs, ex- 



plaining the key puzzle of the very soft phase of the 
complex order parameter. Furthermore, a striking mass 
divergence is predicted that dictates the occurrence of 
quantum critical point observed at ~ 5.2% doping. Our 
model describes most properties of the superfluid very 
well compared with the experiments, and provides new 
insights into several current puzzles. This simple and 
solvable model offers a useful paradigm answering the 
long-standing challenge of describing the underdoped 
cuprates at low temperature. 

Model in the strong coupling limit 

Our simple model is derived in the following steps. 
First, the effective one-particle kinetics of the low-energy 
doped holes is extracted from the experimental spectral 
functions of the 'normal state' at temperature slightly 
higher than Tc^ the so-called 'pseudo-gap' phase. Next, 
we consider the implications of strong coupling limit and 
define a bosonic description of the paired holes. Finally, 
the pivoting motion of the pairs is derived presenting an 
extended hard-core constraint. 

Figure lb shows the dispersion of the main features 
in the measured spectral functions by ang ular-resolved 
photoemission spectroscopy (ARPES)^^^*^, and our the- 
oretical fit via a simple tight-binding Hamiltonian. One 
notices immediately the well-known observation that the 
dispersion is not rigid against the doping level. The 
close resemblance to the previously published t-J model 
solutions^^ in Fig. Ic suggests that this strong doping- 
dependent renormalization originates primarily from the 
competition between the bare kinetic energy and the anti- 
ferromagnetic (AF) interaction. As a consequence, the 
resulting fully-dressed first, second and third neighboring 
hopping parameters, t, t' and t" ^ of the holes (Fig. Id) 
demonstrate strong doping dependence as well. Inter- 
estingly, as 5 decreases, t" is found to increase steadily 
approaching the value of t' ^ and then rapidly exceeds t' 
right at Sqcp! This is apparently not a coincidence, and 
reveals an important clue to the nature of QCP to be 
discussed below. As a reference. Fig. le also shows the 
effective mass of the doped holes, tti^, in three directions. 
In agreement with the experimental observation^, ?72^ 
(red line) is about constant for 6 > 5.2%. 

Note that the resulting hopping parameters should be 
properly understood as only convenient representations 
of the average effective kinetics of one-particle propagator 
under the full renormalization of the many-body interac- 
tions. They certainly do not contain explicit information 
of two-particle interaction nor any decoherent processes 
that broaden the spectrum. Since the fully dressed t is 
negligibly small, understandable from the strong AF cor- 
relation, it will be dropped from our further analysis. 

Now, taking the observation of strong diamagnetism at 
T > Tc as a direct evidence of existence of locally coher- 
ent pairs, let's consider the strong coupling limit, where 
the pairing strength between doped holes is stronger than 
the renormalized kinetic energy, such that at low tem- 
perature T < Tc the doped holes are mostly paired with 
other doped holes. This non-conventional condition is 



well realized in the underdoped region, as made explicit 
from the well-established Gutzwiller projectiorP^ with 
strong AF correlation, and/or with the formation of bi- 
polaron^^. Furthermore, since it is unlikely that doped 
holes can doubly occupy the same site in a weakly doped 
AF Mott insulator, the doped holes are then expected 
to be mostly nearest neighbor to each other under such 
a strong coupling. It is thus convenient to employ a 
bosonic picture consisting of pairs of tightly-bound first- 
neighboring doped holes, 6^- = cLcL located at site i 
and j with opposite spin. In essence, the low-energy 
physics of the underdoped system is then mapped to a 
dilute bosonic gas, a well understood problem that can 
be solved quite accurately. 

The strong coupling limit is the exact opposite to the 
standard BCS theory whose validity relies on the weak- 
ness of pairing interaction in comparison to the kinetic 
energy. Not surprisingly, the key low-energy physics is 
quite different in this case. For example, consistent with 
Anderson's "mammoth vs. mouse" argument^^, typical 
concerns about the pairing mechanism (i.e. exchange of 
glue particles like the phonon or magnon) is irrelevant at 
T <Tc^ since pairing is so much stronger than Tc that the 
high-energy depairing process is of no direct relevance to 
the properties of the superfluid itself. Instead, the key 
factor is the degree of phase coherence, the focus of this 
study. Interestingly, the removal of interaction strength 
from the low-energy physics leaves the effective kinetic 
energy the only relevant energy scale, allowing a simple 
and generic description shown below. 

As another example, the superconducting pairs are 
typically considered as pairs of electrons (or holes) made 
of quasi-particles with opposite momentum, but here a 
tightly bound pair consists of superposition of the whole 
momentum space. In the strong coupling limit, holes be- 
ing part of a pair cannot be considered as quasi-particles 
at all, due to the strong pairing interaction. As shown 
elsewhere^ ^- the very low-energy quasi-particles observed 
by ARPE^SEoMIl ^^ STM'^^ are in fact not the typical 
Bogoliubov quasi-particles that get knocked out of the 
pairs. Similarly, the typical decay of superconducting 
pairs via depairing into two low-energy quasi-particles^ 
is negligible at low temperature here^^. All these seem- 
ingly ant i- intuitive physics only reflects the well-known 
fact that a strongly interacting system may not be simply 
considered as a collection of weakly interacting quasi- 
particles. Extreme caution should thus be exercised in 
such limit before applying existing knowledge and intu- 
ition built on the BCS theory. 

Most importantly, the strong binding between the 
paired holes modifies significantly the kinetics of the bo- 
son, as demonstrated in Fig. 2. Consider a single pair 
(blue and red filled diamonds) located in the fermion lat- 
tice in Fig. 2a. The hopping of each hole is now con- 
strained by the requirement of staying paired. As a re- 
sult, only three potential destinations (empty diamonds) 
for each hole are allowed, two via second neighbor hop- 
ping, t', one via third neighbor hopping, t'^ (Higher 



order processes involving temporary pair breaking are 
negligible in the strong coupling limit.) Converted to 
the bosonic description in Fig. 2b, one finds two inequiv- 
alent bosonic sites in a square lattice, each can hop to 
four first neighboring sites with t\ but only two second 
neighboring sites with t^^ along their designated direction 
(denoted by the direction of the ellipsoids). This incom- 
plete lattice forms the building block or our model. The 
resulting bosonic Hamiltonian reads: 



H = 2_\'^U'^^i' + h.c.[-\-constraint] 



(1) 



where / denotes the bosonic lattice sites and tu^ has value 
t' for the first neighbors and t^^ for the designated second 
neighbors as illustrated in Fig. 2b. 

It is essential to note that the bosonic operators, 6's are 
under a serious 'extended hardcore constraint': b]-b],-, = 
if i = i' OT j = j\ inherited from the Pauli exclusion 
principle of the original fermion operators cl^cl^ = 0. In 
the bosonic lattice shown in Fig. 2b, this constraint for- 
bids occupation of the six potential hopping destinations 
of each boson by another boson. In other words, the 
real-space pairs are correlated via a "polite" policy that 
no one is to get into other's immediate hopping path. 
This constraint can be considered as an infinite short- 
range repulsion that dominates the inter-bosonic inter- 
actions in most cases, and is responsible for stabilizing 
our bosonic system against potential phase separation^^. 
Also note that our model can actually be derived equiv- 
alently through the use of full one-particle propagator 
diagrammatically, but the derivation presented here is 
more intuitive and reveals more insights in the following 
analysis. 

Soft phase and mass divergence of d-wave pairs 

Following the standard approach to dilute bosonic sys- 
tems, let's first examine the solution of the kinetic part 
of the Hamiltonian, as it turns out to be able to answer 
all the main questions of this study already. For conve- 
nience, a unit cell is chosen to include four boson sites 
as shown in Fig. 2b. This choice explicitly allows one s-, 
two p-, and one d-wave superposition within the unit cell, 
and equates the doping level per unit cell of boson to that 
of the standard fermion lattice. Fig. 3a illustrates the re- 
sulting bosonic band structure (band energy Ck vs. crys- 
tal moment k) with doping level {S = 15%) above 5qcp- 
Interestingly, the lowest energy band is composed of d- 
wave symmetry (red color) with a single minimum, where 
Bose-Einstein condensate (BEC) would take place at low 
temperature. One thus expects the superfluid behavior 
(boson flowing without dissipation) once the constraint 
(interaction) is enabled. That is, the d-wave supercon- 
ductivity would appear in the form of bosonic superflu- 
idity in the presence of a true BEC. 

The real-space point of view provides further insights 
into the low-energy Hilbert space of the solution. As 
shown in Fig. 2a, the local Wannier function, calculated 
as Fourier transform of the Bloch functions of the lowest 
band, has clear d-wave symmetry with nodes along the 



(tt, tt) directions of the standard Fermion l attice, in agree- 
ment with the experimental observations J -'^^^^^ J^ (Recall 
that our bosonic lattice is 45-degree rotated from the 
fermion lattice.) That is, the d-wave symmetry is en- 
tirely a local phenomenon unrelated to the phase coher- 
ence. Therefore, the unconventional d-wave symmetry 
only serves as a "form factor" of experiments that ob- 
serve the pairs via their coupling to the fermions. It is, 
however, completely irrelevant to the general behavior of 
superfluid itself. On the other hand, as far as phase co- 
herence's concern (how bosons with the "shape" of the 
Wannier function aligning their phases across the sys- 
tem), the zero momentum of our BEC is of no difference 
from that of the standard BCS s-wave superconductivity. 

Our resulting local d-wave symmetry is entirely driven 
by the screened kinetic energy, consistent with a previous 
study of the t-J model.!^ Indeed, the positive sign of t' 
of the local pair prefers energetically opposite sign of the 
wave function across first neighbors, favoring a (i-wave 
symmetry. On the other hand, the positive sign of t" 
favors opposite sign across the second neighbors, thus p- 
wave symmetry (denoted via green color). This produces 
dominant incoherent local p-wave pairs at 5 < Sqcp^ as 
shown in Fig. 3c. 

The competition between (i-wave and p-wave leads to 
a crucial physical effect, namely a significant mass en- 
hancement. Even near the optimal doping {6 ^ 15%), 
the comparable value of t" and t' leads to a large ef- 
fective mass m* = {h? /P)d^ek/dk'^ ^ 12m^ ^ 59me 
(with lattice constant /), corresponding to a penetration 

depth A = A / 4^1? n ^ 7000 A, in reasonable agreement 

with the experimental valu^^ (taking Uq ~ 5/2 per Cu 
atom). Furthermore, as 5 approaches Sqcp^ t" grows 
to the value of t' , reducing the separation of the d-band 
and the p-band, and in turn flattening the (i-band. The 
effective mass of the d band thus increases significantly 
(Fig. 3d), and eventually diverges at the QCP, where p- 
and d-wave become degenerate at t" = t' (Fig. 3c). Since 
the superfluid stiffness correlates with the inverse of the 
effective mass, the large mass enhancement would dictate 
a very soft phase. 

Our results also lead to an unexpected exotic conclu- 
sion: The QCP at the end of the underdoped superconduc- 
tivity region is dictated by a mass divergence of the local 
pairs. With an infinite mass, the phase of local pairs be- 
comes entirely free even within short distance, since they 
can no longer communicate with each other. Even at 
zero temperature, this complete lost of phase coherence 
would suppress entirely the BEC. Equivalently from the 
bosonic perspective, the infinite number of minimum in 
the flat band dispersion in Fig. 3b can no longer support 
a BEC. Consequently, the diamagnetic response would 
also be strongly suppressed near the QCP due to lack 
of short-range coherence, in perfect agreement with the 
experimental observation^^. 

Interestingly, at J < 5qcp-, the mass of the p-wave 
pair still diverges along the anti-nodal (tt, tt) directions. 



dictated by the symmetry of our lattice in Fig. 2b, ren- 
dering BEC impossible as well. It would be curious to 
investigate the scattering against entirely incoherent p- 
wave, instead of the impurities, as a microscopic origin 
of the well-known glassy behaviors found in this region. 

Note that in great contrast to the mass divergence as- 
sociated with the metal-insulator transition^^, the above 
divergence results from the competition between kinetic 
effects of t' and t" . This competition can be considered 
as a destructive interference resulting from the wave na- 
ture of the quantum particles. As illustrated in Fig. 4, in 
comparison to the weakly bound superconducting pair in 
the BCS theory, a tightly bound local pair is subject to a 
strong interference effect, and thus reduced net kinetics 
(or enhanced mass). The surprising mass divergence at 
the QCP can thus be viewed as a complete destructive 
interference that renders the hole pairs immobile. 

Our predicted mass divergence near QCP is strongly 
supported by indirect experimental evidences. The in- 
verse square of the experimental penetration depttP^, 
A~^, of YBa2Cu302y (YBCO) shows non-linear behav- 
ior against doping (Fig. 5a). More significantly, near the 
QCP, measurement performed on a single YBCO sample 
under repeated annealing to modulate the doping leveP^ 
cleanly shows a strongly non-linear correlation between 
A~^ and Tc (Fig. 5b). The zero slope at A ^ in both 
cases can be considered indirect evidence of the mass 
divergence. Nonetheless, a direct measurement of the 
superfluid inertia would be highly valuable to unambigu- 
ously verify our conclusion. 

In-plane superfluid properties 

So far, our model has managed to address the key is- 
sues of this study without use of any free parameter. Be- 
low we verify the applicability of our model in describ- 
ing the in-plane behavior of the superfluid in compari- 
son with the available experimental data. To this end, 
we adapt numerically the well-established self-consistent 
T-matrix approach^^ ^^ to the screened interaction and 
apply the Bogoliubov transformation^^ to study the su- 
perfluid properties (see supplementary materials for de- 
tail). Taking the extended hardcore constraint as the 
leading "interaction" corresponding to a hard-sphere ra- 
dius comparable to the lattice constant /, the scattering 
length of the low-energy particle (the Wannier function) 
is estimated to be a ~ 1/4. 

Naturally, a true BEC also requires tunneling of the 
bosons out of the plane. This is, however, not easily 
treated accurately due to the highly "effectively disor- 
dered" nature of the strongly coupled pairs in our pic- 
ture. Given the weak inter layer hopping of the holes 
{tz ~ 12meV^^), the out-of-plane tunneling of a pair may 
encounter one of the following three possibilities. First, 
one of the holes in the pair hops to the next plane and 
happens to encounter another unpaired hole to form a 
new pair. This first order process, however, is unlikely 
due to the rare supply of unpaired holes in the strong 
coupling limit. Second, one of the holes hops to the next 
plane, followed by its partner. Taking the magnetic cou- 



pling J/4 ^ 33me]M^ as the representitive pairing inter- 
action [/, this second-order process gives tunneling coeffi- 
cient tj_ = 2tl/U ~ 9meV (or smaller if other real-space 
attractions also contribute). Of course this second pos- 
sibility is only allowed when the spins across the plane 
happens to align ferromagnetically in the vicinity of the 
pair. This brings the third possibility that the tunneling 
would be forbidden by the AF arrangement across the 
plane. Obviously, such a violent (and strongly ^-, T- and 
site-dependent) on/off modulation of the tunneling is ex- 
pected to limit significantly the superfluid density in the 
out-of-plane direction, through shorter bosonic lifetime 
or possibly development of localization^. This analy- 
sis also indicates that while the in-plane AF correlation 
might have been one of the leading factors to provide a 
strong pairing of holes, a strong AF correlation across 
the plan is actually very damaging to the formation of 
BEC and to the out-of-plane superfluid density. 

We thus limit the following discussion of the superfluid 
to the in-plane properties only. For simplicity, a fixed 
bosonic tunneling tj_ ~ 3meV is chosen as a parameter 
that resembles a rough average over the above three pos- 
sibilities. This should be understood as a crude estimate 
to reproduce the average degree of coherence only in the 
study of in-plane superfluid density. (Such a simple aver- 
age is certainly inadequate to give anything sensible for 
the out-of-plane superfluid density.) All the results re- 
ported below corresponds to this single set of parameters 
t± and a. At low temperature, the scattering with the 
BEC results in quasi-particle energy Ek = \/el + 276, 
consisting of linear dispersion (sound propagation) at 
the long wave-length for energy within the characteris- 
tic scale 7. Therefore, the flow of the pairs suf fers no 
dissipation below the superfluid velocity, Vs = h^J^|vfl^ . 

Figure 5c shows the resulting phase diagram, in which 
Tq is determined by solving J DOS{w) ^ fB{w)dw = n 
numerically, where n = S/Qq is the hole carrier density 
and DOS{w) is the "bare" density of state correspond- 
ing to the band structure in Fig. 3. Owing to the above 
choice of parameters, the resulting Tq reproduces the ex- 
perimental phase boundary quite well at the underdoped 
region. (The experimental Tc is expected to deviate from 
our Tq near the optimal doping, where amplitude fluctu- 
ation becomes relevant.) 

Moreover, all our key conclusions reflects explicitly in 
the superfluid density, n^ = n — n^, where the normal 
component of the number density n^ is obtained via n^ = 
lim^^^ = jDOE\w)fB{w)dw where DOE{w) ^ 

/ 0S"^ (^ - ^k) = {2we/3{e + j))DOS{e) with e = 
Y^7^ + w'^ — 7. Indeed, the zero-temperature superfluid 
density in Fig. 5d vanishes rapidly near the QCP, and re- 
mains zero below QCP, despite the finite density of local 
pairs (red dotted line.) Interestingly, our resulting pene- 
tration depth compares quite well with the experimental 
trends against S and Tc (Fig. 5a and 5b.) 

Several important issues concerning the temperature 
dependence of superfluid density can now be addressed 



by our results. As shown in Fig. 5e, experimental ns{T) 
at ^ ^ 7% shows a approximately linear reduction at 
low temperature, in great contrast to the exponentially 
small reduction in standard s-wave BCS theory. This 
experimental observation was originally viewed as an ev- 
idence of the d-wave structure of the superconducting 
order, ^^ resulting from superfluid decaying into quasi- 
particles near the nodes, where the superconducting gap 
decreases to zero. Also shown in Fig. 5e, such rapid re- 
duction of ris (T) can be reproduced very well by our cal- 
culation. As pointed out earlier, in the strong coupling 
limit assumed in this work, decay into quasi-particles is 
inefficient. Instead, in spirit similar to the alternative 
explanation by Emery and Kivelson,!! the rapid reduc- 
tion of ns{w) in our results originates from the thermal 
depletion of the BEC, completely unrelated to the local 
d-wave structure. 

A more serious question was raised in the recent mea- 
surement of A near the QCP, also shown in Fig. 5e, 
namely the absence of presumable quasi-2D nature of 
the superfluid. The low-energy hol e carr iers in cuprates 
are well accepted to be quasi-2D.^^^^ Naturally, the 
sharp drop of ns{T) at T - O.ST^ for S ^ 7% was 
commonly interpreted as the signature of the so-called 
BKT transitioiP^H^ associated with the quasi-2D na- 
ture. However, such rapid reduction is entirely absent 
near QCP. In addition, in both 7% and near QCP, ns{T) 
appears to be linear near Tc, casting doubt on the quasi- 
2D assumption of the superfluid. 

Our results again capture these features very well and 
thus provide a potential resolution to this important is- 
sue. It turns out that the sharp drop of ns{T) at higher 
doping (e.g. 15%) only reflects approximately the energy 
scale of 27(T), the energy scale of linear dispersion of 
the bosonic quasi-particles. Throughout the underdoped 
region, the enormous in-plane mass enhancement sup- 
presses significantly the quasi-2D nature of the superfluid. 
Especially close to the QCP, the near divergence of the 
in-plane mass obviously will not support any quasi-2D 
behavior, even though there still exists qualitative dif- 
ference in the out-of-plane direction due to its strongly 
disordered characteristics. 

Finally, our model gains further verification from com- 



paring experimental critical current with the intrinsic 
theoretical limit jc <= eUsVs. Figure 5f shows three 
measurements^S^*^ all of which consists of an reflection 
point around T ^ 0.8Tc, also produced nicely in our cal- 
culation. 

As final remarks, our main findings, namely 1) kinetic 
origin of the (i-wave symmetry, 2) soft phase due to mass 
enhancement and 3) QCP dictated by mass divergence, 
all depend only on the ARPES measurement and require 
no free parameter. Our current study does not include 
the special case of the 'stripe' phase^^, although it should 
be possible to describe it as a long-range order of the lo- 
cal pairs, with a microscopic derivation of the anti-phase 
'd-density wave' phenomenon. ^^^^ Our bosonic picture is 
built on bond-centered boson and is thus ideal for the de- 
scription of the observed broken C4 symmetry in STM^^. 
Finally, extension of current work to include pairing be- 
yond first neighbors (beyond the strong coupling 'limit') 
can be formulated in a similar manner, but such exten- 
sion will not change the symmetry properties of the solu- 
tion, and thus should preserve qualitatively the key find- 
ings of this work. 

Conclusion 

In summary, utilizing the strong AF correlation, we 
propose a very simple microscopic model of strong 
coupling limit of superconductivity in the underdoped 
cuprates, where physics is evidently dominated by phase 
fluctuation. The d-wave symmetry is found to be en- 
tirely a local phenomenon driven by the kinetic energy, 
and the observed superconductivity can be understood 
as a superfluid of a dilute bosonic gas of local pairs. Fur- 
thermore, a severe competition between ffist and second 
neighbor hoppings leads to a colossal mass enhancement 
at underdoped region, which in turn produces the ob- 
served very soft phase. The presence of quantum critical 
point at the end of the underdoped superconducting re- 
gion is predicted to be dictated by a mass divergence of 
the pairs with which the phase coherence is totally lost. 
Our results answer several key questions on the under- 
doped cuprates and appear quite adequate for describing 
the in-plane superconductivity in the underdoped region 
in a simple and coherent manner. 
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FIG. 1: Fig. 1 (a) Phase diagram of LSCO in temperature 
T vs. hole doping (^^). Tmf gives the Tc from mean-field 
theor}/^^ . Te denotes the temperature of appearance of phase 
coherence, and Tc the experimental superconducting transi- 
tion temperature, (b) Illustration of our fitting to the exper- 
imental band dispersion, (c) Band structure of t-J modeF^. 
(d) Parameters t, t' and t" of holes obtained from (b). (e) 
Doping dependence of carrier mass, ?7i^, in different direc- 
tions indicated by the arrows in the inset. 



FIG. 2: Fig. 2 (a) Illustration of kinetic hopping of a pair of 
holes (filled squres). Open squres indicate the allowed des- 
tinations under the pairing constraint. Ellipsoids denote the 
same processes from the perspective of the pair as a whole. 
(b) Illustration of bosonic kinetic processes of the paired holes. 
The solid/dashed lines denote first- /second-neighbor hopping 
paths of strength t' /t" . The yellow area shows the region of 
the 'extended hardcore contraint' where occupation of aother 
pair is forbidden. The dotted square gives the unit cell for 
our calculations. 



FIG. 3: Fig. 3 The band dispersion of kinetic-only solution at 
5 > ^QCp(a), 8 = ^QCp(b), and 5 < Sqcp{c). Red, green and 
blue color represents the local d-, p-, and s-wave symmetry. 
Insets in (a) and (c) illustrates the symmetry of the Wannier 
function corresponding to the lowest set of bands, and the 
kinetic effects that drive the symmetry, (d) Effective mass of 
the pairs, m*, and the measured holes, m^. 



FIG. 4: Fig. 4 Illustration of heavier mass resulting from 
stronger binding. 



FIG. 5: Fig. 5 (a) Non-linear doping dependence of penetra- 
tion depth, (b) Non-linear temperature dependence of pene- 
tration depth, (c) Phase diagram with our model (blue dots; 
other symbols identical to Fig. la), with inset showing detail 
near the QCP. (d) Doping dependence of the superfluid den- 
sity, (e) Temperature dependence of superfluid density for 
various doping levels, (f) Temperature dependence of critical 
current. 
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